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Abstract 

The symmetry analysis of Ermakov systems is extended to the gen¬ 
eralized case where the frequency depends on the dynamical variables 
besides time. In this extended framework, a whole class of nonlinearly 
coupled oscillators are viewed as Hamiltonian Ermakov system and 
exactly solved in closed form. 


1 Introduction 

The first study of an Ermakov system was published by the end of the last 
century |l| , but the motivation for a detailed investigation of Ermakov sys¬ 
tems intrinsic properties, their proper characterization and their application 
scope derived mainly from the works of Lewis who rediscovered them in 
the sixties [|], and of Lewis and Riesenfeld, who devised one of their most 
important application [Q, the quantization of the time-dependent harmonic 
oscilator. 

In its original formnlation, the Ermakov system was nothing bnt a time- 
dependent harmonic oscillator conpled to Finney’s 0 eqnation. Later on, 
Ray and Reid extended the basic concept by incorporating two extra 
arbitrary fnnctions in the original equations. Snch generalization highly im¬ 
proved the reach of the Ermakov systems with no harm to their key property, 
namely the existence of a constant of motion, the Ermakov-Lewis invariant. 
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Early developments dealt with Ermakov systems where the frequency [Q 
depended solely on time. In a subsequent generalization 0, the dynami¬ 
cal variables were allowed to participate in the argument of the frequency. 
The resulting generalized systems have rather interesting properties and po¬ 
tential applications ii, but has stayed mostly as a fertile field for further 
prospects. In particular, the Lie point symmetry analysis of Ermakov systems 
performed by Leach and co-workers @-|12| considered frequency functions 


depending at most on time. Inspired by the results of Leach and co-workers, 
we reconsider the issue of geometric symmetries for generalized Ermakov sys¬ 
tems where the frequency is a function of both the dynamical variables and 
their first derivatives besides time. By doing this we hnd a much wider 
class of dynamical systems that aside from possessing the Ermakov-Lewis 
invariant admits a Lie symmetry in a direct and natural generalization of 
that proposed by Leach. As a first usefull application of this extended view¬ 
point we show that some systems which would otherwise not fit in the older 
framework are readily treated as Ermakov systems. 

This letter is organized as follows. In section II, the class of generalized 
Ermakov systems with Lie point symmetry is determined. In Section III 
a class of coupled nonlinear oscillators are shown to constitute a general¬ 
ized Ermakov system that possesses geometric symmetries and are exactly 
solvable in closed form. The conclusions are presented in section V. 


2 The Symmetry Analysis 

All Ermakov systems, both usual and generalized, can be represented Q] by 
the pair of equations 

x + n'^x = -^F{y/x), ( 1 ) 

yx^ 

y + n^y = 0 , ( 2 ) 

where F is an arbitrary function of the indicated argument and hi, hereafter 
called the generalized frequency, is an arbitrary function of the time, the 
coordinates and their derivatives. For physical reasons, we restrict the study 
to the cases where the frequency depends at most on the velocities, that is, 
we take hi = r2(t, x, y, i, y) only. 


2 





The pair of equations TO encompasses the whole class of Ermakov 
systems in two spatial dimensions. In the traditional notation, the Ermakov 
systems were written in terms of three arbitrary funtions and not only 
two as we have proposed [§ elsewhere. This more concise representation is 
possible thanks to the generalized character of The usual representation 
is obtained whenever is chosen in the form 

- ^g{x/y), (3) 


where g{x/y) and to are arbitrary functions of the indicated arguments. 

The use of generalized frequencies considerably expands the reach of the 
Ermakov system concept. Among the systems that can be treated as partic¬ 
ular Ermakov systems in the generalized form, we quote the Kepler-Ermakov 
system m and the Lutzky’s integrable system ||14|| . 


In order to perform the Lie symmetry analysis of generalized Ermakov 
system m we follow the approach of Govinder and Leach for the usual 
case BUI, and consider initially the symmetries of the equivalent system 


xy-yx + —F{y/x) = 0 , 
x^ 

y + n‘^{yi,±,t)y = 0, 


(4) 

(5) 


where the notation x = (x, y) was introduced. The advantage of dealing first 
with the equivalent equation (|) derives from its independence of G. This 
feature provides a straightforward way of finding its symmetry generator, 
which is basically that obtained by Govinder and Leach in the symmetry 
analysis of the usual Ermakov systems. 


G = y|+ppx. 


A 


( 6 ) 


12 


where p{t) is an arbitrary differentiable function. As pointed out in ref. 
the Ermakov first integral I is invariant under the first extension of G. 

For the complete Ermakov system, that is, for equations (M to obey the 
same symmetry, some restriction must be imposed on the functional form of 
G. This restriction is determined by requiring the invariance of Eq. (|^) under 
the action of the first extension of G. We apply this criterion, and obtain 
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the most general frequency that preserves the invariance of the generalized 
Ermakov system under the action of G, 

=-^ + \a{x/p,pk-p^), ( 7 ) 

p p^ 

where a is an arbitrary function of the indicated arguments. Notice that 
velocity dependence is also possible, a fact that considerably expands the 
class of Ermakov systems invariant under the symmetry transformations of 

G. 

To summarize the results so far we note that with no a priori assump¬ 
tion on the functional form of the frequency, we hnd that the most general 
Ermakov system possessing Lie point symmetry must be of the form 


X + I 

( cr 

1 d 

-- 

1 X = -\F{y/x) , 

(8) 



P) 

yx^ 


y + 

4 

_ P^ 

\y = 0. 

(9) 


\p 

P) 

F 


where a is arbitrary but restricted 

to 



a 

= ct(: 

t/p, p±- p^). 

(10) 


The associated transformation group is that generated by G. 

Let us show now that the whole class of usual Ermakov systems are 
recovered as special cases of the present treatment. For this purpose we use 
the relation 

p + u,'\t)p=^, ( 11 ) 

to define a new function uj{t) in terms of p and the constant, but otherwise 
arbitrary, flo- Despite its appearance of a differential equation, (0) can be 
alternatively considered as the dehnition of uj in terms of p. Let us adopt 
this viewpoint and also dehne a transformed a 

d(x/p, px - px) =-Do + — -g{xly) + a, (12) 

xy^ 

where we introduced, for convenience, a spurious but otherwise innocuous 
function g{x/y), which can be absorbed in the dehnition of a. These redeh- 
nitions transform the generalized form of the Ermakov system (^|-^ into its 
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more traditional form 


X + + ^)x = 

y + = 


In equations (|I3| -|T^) we replaced 


F{y/x) = fiy/x) + 


-^f{y/x) , 

yx^ 

(13) 

F,<Axly). 

(14) 


(15) 


in order to facilitate the comparison with previous results on the Lie group 
of usual Ermakov systems. We should stress, however, that the superfluous 
function f{y/x) was introduced only in order to reconstruct the traditional 
form. 

We now notice that the pair of equations ([T3|-[I^ reduces to the usual 
form with simple time dependence in the frequency only when ct = 0. For 
this particular class of frequencies the quasi-invariance transformation 


X = x/C y = y/C J 


removes oo 


16| from the equations of motion, which become 
1 


x” = 


;f{y/x ), 


y" = ^9{x/y ), 


(16) 


(17) 


yx^ xy^ 

whenever C{t) is a solution of the time-dependent harmonic oscillator 

C + uj^{t)C = Q. (18) 


The symmetry group of the transformed Ermakov system (1^) is SL{2,R) 
as already shown in 0 . This however is not the case with equations (|T3|-^) 
which are more general than the Ermakov system with SL{2, R) symmetry. 
The extra symmetry freedom of the system (|T^|^) derives from a, which can 
depend on the coordinates or the velocities. Notice however that the trans¬ 
formations belonging to the SL{2, R) group are among the transformations 
generated by G. 

The symmetries considered so far do not require any particular form for 
the function a, nor depend on any preliminary coordinate transformation 


5 








like (|^). In this respect, it should be noted that, despite its apparent sim¬ 


plicity, the transformation (^) is not well defined when a; depends on the 
dynamical variables, a situation in which equation (|^ becomes meaningless. 
In addition, even for uj = the time-dependent harmonic oscillator (|1^) 
does not always possess a closed form solution, a fact that may preclude its 
application in many interesting practical situation. Similar remarks apply 
to the linearization of the Ermakov system proposed in [Q, which starts 
exactly with the transformation (^). Needless to say, such a linearization 
process applies only to the simplihed system ([T^ and not to the generalized 
Ermakov system (||-^. However, it should be stressed that other generalized 
Ermakov systems involving extensions to multi-components and higher di¬ 
mensions have recently been introduced which do have the underlying linear 


structure of the Ray-Reid system [17, IS 


3 Hamiltonian Generalized Ermakov Systems 

The Ermakov systems with SL{2, R) symmetry in two spatial dimensions 
are solvable in closed form |^. The more general Ermakov systems 
however, are not always exactly solvable and this, certainly, is the price paid 
for the extra generality. Some additional structure, perhaps a Hamiltonian 
formalism, must exist in order to reobtain equations that allow for exact 
solutions. To present one application example and to gain insight in this 
respect, let us consider the two dimensional oscillator with nonlinear coupling 


introduced by Ray and Reid [19 


x + u‘^{t)x = xp ^G{xy/p^), 
y + uj'^{t)y = yp-^G{xy/p‘^), 

where p{t) is an auxiliary function satisfying Finney’s equation 

p + uj^{t)p= 1/p^. 

Ray and Reid found a first integral for namely 

xy 

J ^ {px - fix){py - py) + — - J GirjdT, 


( 19 ) 

( 20 ) 


( 21 ) 


( 22 ) 


6 















but did not obtain its complete solution. In this section we show that the 
generalized Ermakov system approach combined with the symmetry analysis, 
yields the corresponding complete solution, in closed form. 

Let us hrst show that (p!9| -pOD is in fact a generalized Ermakov system. 
This is easily verihed by choosing in (@-^, F = 0 and a function of the 
coordinates 

- -^Gixy/p-^). (23) 

From this viewpoint the amplitude of the angular momentum 

l = \ixy-yxf, (24) 


becomes the corresponding Ermakov-Lewis invariant. Moreover, 
form a Hamiltonian Ermakov system, a property that has already played 
a major role in several applications 
Hamiltonian is 


20 


For the present case the 


rxy/p'^ 


H = PxPy + uj'^xy - / G{t) dr , 

P 


(26) 


which is not a constant of motion, since the system is non-conservative. 

The frequency (PB|) is of the right form ( 0 ) and therefore guaranties the 
invariance of ([l9|-pOD. The use of the auxiliary equation to substitute p 
for transforms into 


iV‘ = -'- + -^(l-G(xylp-‘)), (26) 

P P^ 

which is indeed of the general form (^. 

The Hamiltonian structure and the symmetry can now be further ex¬ 
ploited to reduce the problem to quadrature. For this purpose, the use of 
canonical group coordinates. 


u = x/p, v = ylp, T 


p dt, 


(27) 


is an essential artihce that transforms the symmetry generator into a simple 
time translation, G = d/dT. In this new variables, the equations of motion 
become the autonomous system 


u" + u = uG{uv ), v" + v = vG{uv ), (28) 


7 



where primes denote derivatives with respect to the new time. Moreover, 
this system is described by the Hamiltonian 

/ UV 

G{t) dr, (29) 

which is the transformed version of the hrst integral J given by The 

Hamiltonian K in the canonical group coordinates is quadratic in momentum. 
In such cases p, the appropriated variables for quadrature are 

q = UV , s = v/u . (30) 

In these coordinates, the constants of motion become 

= qs'/s, (31) 

J = (g'^ - 2/)/4g + g - y G{t) dr . (32) 

The quadrature of the last two equations give successively g(T) and s{T). 
The map back to the original variables is a pure algebraic task of using 

x‘^ = p^q/s, y^ = p‘^qs, (33) 

to substitute {x,y) for (g, s). One hnal step is the determination of T(t) 
which can be obtained by the integration of the last equation in (pTl) . In 
the quadrature sense, this completes the integration process. Recall that 
all the necessary structure and the correct choice of variables was dictated 
by the Hamiltonian character, the generalized Ermakov form of the system 
and by the presence of symmetry. In particular, the analytical form of the 
solution can be found explicitly in terms of elliptic functions when G{t) = 
Ci/r^ + C 2 + CsT + C 4 r^, where c* are arbitrary constants. 


4 Conclusion 

Generalized Ermakov systems (^-^ are of a very general nature, and their 
complete application scope is still not fully outlined. In this paper we deter¬ 
mined their Lie point symmetry group and showed that a set of nonlinear 
coupled oscillators belonging to this class of systems is in fact explicitly solu¬ 
ble in term of quadratures. This provides a stimulating preview of the variety 


of applications in which generalized Ermakov systems may play an important 
role. 

To conclnde, we mention that Lie symmetries of the type generated by 
and some modifications of it occur in several problems in Physics. Some 
non Hamiltonian systems with velocity dependent frequency, a whole class of 
two and three-dimensional charged particle motion, some perturbed magnetic 
monopole and a time-dependent Kepler problem fall in the category of such 
systems and are frequently amenable to exact solution in closed form. 
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